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Owing to their fractional electric charges, quarks are blind to transformations that combine a 
color center phase with an appropriate electromagnetic one. Such transformations are part of a 
global Zg-like center symmetry of the Standard Model that is lost when quantum chromodynamics 
(QCD) is treated as an isolated theory. This symmetry and the corresponding topological defects 
may be relevant to non-perturbative phenomena such as quark confinement, much like center 
symmetry and ordinary center vortices are in pure SU(A^) gauge theories. Here we report on 
our investigations of an analogous symmetry in a 2-color model with dynamical Wilson quarks 
carrying half-integer electric charge. 
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1. Introduction 

In spite of tremendous research efforts, the phase diagram of QCD matter at non-zero temper- 
ature and chemical potential is not well understood. There remains much to learn about the exact 
nature of the transition of hadronic matter into a deconfined plasma in the hot and dense conditions 
of heavy ion collisions and the early universe. 

For the finite temperature deconfinement transition of infinitely heavy quarks, at least, the 
importance of center symmetry has been established [1]. The deconfinement transition of pure 
SU(A'^) gauge theory in (i+ 1 dimensions mirrors the Zj^2> symmetry breaking of a (i-dimensional Zj^ 
spin model [2], with spacelike center vortices playing the role of spin interfaces. Their percolation 
at low temperatures leads to confinement by disordering the Polyakov loop, while their dynamical 
suppression at high temperatures allows for an ordered, deconfined phase [3, 4]. The connection 
with spin models is especially striking for the second order phase transitions of SU(2) and SU(3) 
in 2-1-1 dimensions [5, 6]. Here the self-duality of the corresponding 2-dimensional spin models 
is manifested in the universal behavior of center vortex and electric flux free energies, which are 
reflections of one another about criticality [7]. 

The role of center symmetry is, however, obscured by the inclusion of dynamical quarks in 
the fundamental representation. Center symmetry is explicitly broken and the finite temperature 
transition of QCD becomes a smooth crossover at zero chemical potential [8, 9]. 

This neglects the electric charge of quarks, which is commonly expected to only require small 
perturbative corrections. Note, though, that the inclusion of the quarks' fractional electric charge 
leads to a global Ze-symmetry that combines the centers of the color and electroweak gauge groups. 
The physical realization of this symmetry could have non-trivial implications for confinement and 
the phase structure of the Standard Model. In these proceedings, we report on our ongoing study 
of such a symmetry in 2-color QCD with half-integer electrically charged quarks. 

2. Hidden symmetry 

We will give only a brief overview of the global center symmetry of the Standard Model and 
its analog in our toy lattice model. The interested reader is directed to our proceedings from last 
year's lattice conference [10], as well as [11] for a review of the group structure and [12] for a 
different lattice realization of the symmetry. 

Since quarks carry fractional electric charges Q = |e or —\e, the color and electromagnetic 
phases in the pair of combined transformations 

^^iK/^^^lTiQle^^ ^^-ilnl^ ^^-ilTtQle^ ^ 3U(3) ^ U(l),„, (2.1) 

cancel precisely. What's more, they act trivially on all other particles in the Standard Model, which 
ai^e blind to the center of SU(3) and caiTy integer electric chage. Electric charge is related to 
hyperchai^ge Y and the third component of weak isospin ?3 by Qje = tj, -\-Y /2, with e'^^^'^ = — 1 G 
SU(2) and Y quantized in units of 1/3. The symmetry is therefore generated by 

(e'2''/^-l,e'''^) G SU(3) xSU(2) xU(lV, (2.2) 

which gives six elements including the identity. 
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A global symmetry brings with it the possibility of a phase transition characterized by spon- 
taneous symmetry breaking. In this case, one expects the transition to be driven by topological 
defects (vortices) that carry both color and electromagnetic flux. That is, center vortices with an 
additional electromagnetic Dirac string. Could this be relevant to (de)confinement? 

3. A toy 2-color world 

Our starting point is 2-color QCD plus electromagnetism with 2 flavors of Wilson fermions in 
3 + 1 dimensions. By including 'up' and 'down' quarks with fractional charge zt^e relative to the 
U(l)em gauge action, we obtain a model with a global Z2 symmetry. The lattice action is 

S=-Y^( ^Re Tr Up + j8„„ cos 9^ + , (3. 1) 

plaq. V ^ / 

where Sf,w is the usual Wilson fermion action with the distinction that parallel transporters for 
quarks are products of an SU(2) color matrix and a U(l)em phase, of the form 

i7^(x)e'^''W/2^ Ui,{x) G SU(2), d^{x) G (-27r,27r]. (3.2) 

The SU(2) plaquettes Up and U(l)em plaquette angles Qp are formed from U^i and 0^ in the usual 
way. In this model, 'fractional charge' means that the parallel transporters for quarks contain half 
the U(l)^m angle relative to the 0^ that appear in the plaquette angle dp. That is, an = —\ 
electromagnetic link for quarks appears as an g'^'' = +1 link in the gauge action. An important 
point is that the 'size' of the compact U(l)em is determined by the quark, which is the smallest 
quantum of electric charge. The range of 6^ (x) is chosen such that we integrate over all possible 
electromagnetic transporters for the quarks, amounting to a double counting for the \]{\)em gauge 
action and integer charged particles. This is consistent with the premise that our compact U(l)em is 
the result of symmetry breaking in a SU(3)— 7>SU(2)xU(l)^m unified theory (see e.g. [13]). 

It is clear from Eq. (3.2) that a combined color center phase — 1 G SU(2) and electromagnetic 
phase e'^i^/^ = — 1 acts identically on the quarks. The model therefore retains a combined color 
and electromagnetic Z2 center symmetry, despite the introduction of dynamical quarks. This is the 
analog of the hidden Zg symmetry of the Standard Model. 

4. Simulations 
4.1 Z2 disorder 

Our preliminary results in [10] indicated that the usual ordering of color links by dynamical 
quarks is negated in our model by the inclusion of fractional electric charge. This is understood by 
analogy with spin systems. In QCD alone, terms in a loop expansion of the fermion determinant 
that wind around the temporal direction favor the center sector in which the traced Polyakov loop 
Pcoi = 1 . They break center symmetry and lead to ordering in much the same way as spins coupled 
to an external magnetic field. The inclusion of fractional electric charge in our model bestows 
quark loops with an additional U(l)em phase which may undo this effect. Consider, for example. 
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the hopping expansion of our Wilson fermion action for Nt = 4 time slices to fourth order in the 
hopping parameter K, 

5/,eff = -16Jc4 ( £cos^-ReTrt/p + 8£cosPe/2-Re/'™/ J + •.., (4.1) 

\plaq. X J 

If the \J(l)em Polyakov loop angle for quarks Pe/2(x) = L|1o' ^t{t,^)/2 is disordered, then the 
Pcoi = 1 sector is no longer favored and SU(2) center symmetry is dynamically restored. Since the 
parallel transporters for quarks possess a Z2 degree of freedom that the \J{l)em gauge action is blind 
to, this is possible even in the Coulomb phase for integer electric charges. The effect is analogous 
to placing an Ising model in a fluctuating magnetic field, or to the Peccei-Quinn mechanism [14], 
in which the coupling of the CP violating term in QCD to an axion field allows for the dynamical 
restoration of CP symmetry. 

Indeed, we found in [10] that the SU(2) Polyakov loop in our toy model was indistinguishable 
from the quenched SU(2) result at values of the hopping parameter that would otherwise cause 
a significant amount of ordering in standard 2-color QCD (k" = 0.15, 0.175). For random, 'hot' 
starts, where the U(l)e,„ links are totally disordered, this was true both for the confined j3e„, < 1.01 
and Coulomb phases /3em ^ 1-01 of the U(l)em gauge action. Since the U(l)e,„ gauge action is bhnd 
to Z2 disorder as seen by quai^k, i.e. e'^f = ±1, it is unable to remove it. 

On the other hand, Z2 disorder was not able to be generated deep in the Coulomb phase for 
'cold', ordered starts (11^ = 1, e'^^/^ = i)_ Here we found the usual behavior for 2-color Wilson 
fermions, which we attibuted to the ergodicity limitations of our Hybrid Monte Cai^lo (HMC) al- 
gorithm. In the Coulomb phase, the many simultaneous fluctuations in the U(l)e,„ angles 
proposed by the HMC update are strongly suppressed. This prohibits the transformation of, e.g., 
an e'^f /2 = i jjnk to an e'^f = — 1 link, even when they are equivalent in the \J{l)e„, gauge action. 

We have since ameliorated this problem by including additional local updates. The first of 
which is a simultaneous transformation of SU(2) and U(l)e„, links, — —Ufi,;e'^^'^^ — —e'^>^^^. 
The only part of the action that is affected by this combined Z2 update is the SU(2) Wilson plaquette 
term. The U(l)e„, gauge action is blind to such phases by default, and the quark determinant is 
also blind because this transformation belongs to our Z2 symmetry. It introduces a thin color- 
electromagnetic vortex around the transformed link. Between HMC trajectories we propose such 
an update on a random link and accept/reject using a Metropolis check for the SU(2) gauge action, 
repeating many times. This is efficient when ^coi is small, but the acceptance rate is exponentially 
suppressed as the average value for the SU(2) plaquette increases with pcoi- It becomes necessary 
to also include updates that flip only a U(l)em link, e'^**/^ — )■ — e'^''/^. The fermion action requires 
an expensive recalculation, but the acceptance rate is much improved for regions of parameter 
where ^coi is large and K is moderate. Given sufficient thermalization, 'hot' and 'cold' starts show 
complete agreement when such local updates are included. This is demonstrated in Fig. 1 , where 
we plot the traced SU(2) Polyakov loop in the U(l)em Coulomb phase at P^m = 2 using the same 
K" = 0.15 as in [10]. It is indistinguishable from the pure gauge result, showing the dynamical 
restoration of SU(2) center symmetry. 
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Figure 1: Volume average of the traced SU(2) 
Polyakov loop on 8^ x 4 lattices with k: = 0.15 and 
Pern = 2. The 'quenched' and 'dynam SU(2)' results 
are from 2-color QCD without electromagnetism. 
'Hot' and 'cold' starts are in the prescence of quarks 
with fractional electric charge. They are consistent 
with the quenched SU(2) result, demonstrating that 
SU(2) center symmetry is dynamically restored. 
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Figure 2: Meson masses on 8^ x 16 at j3„,/ = 1.7 
and jSgm = 2. 'SU(2)' results are from standard 2- 
color QCD. Inclusion of a compact \J{l)coi and frac- 
tional electric charge dramatically shifts the mass 
scale. Much larger values of K are required in order 
to obtain an equivalent ratio of p and n masses at a 
given pcoi- 



4.2 Mass scales 

The values for the hopping parameter K used in our test simulations were based on the 2- 
color results of Skullerud et al. [15]. Given the large impact that fractional electric charge has 
on the SU(2) Polyakov loop, however, it is not reasonable to expect the mass scale to remain un- 
changed and for these K values to remain appropriate. As such, we performed zero temperature 
measurements of meson masses on 8^ x 16 lattices at (C = 0.178, ficoi = 1-7, and jS^m = 2 to di- 
rectly compare with [15]. Masses were extracted from all-to-all 2-point coiTclation functions with 
pointlike sources, using dilution in the spin and time indices in accordance with [16]. 

The inclusion of fractional electric charge has a dramatic effect, even though we are in the 
Coulomb phase for the U(l)e„, gauge action. At these parameters, the 7i and p masses in lattice 
units are roughly twice those of the standard 2-color result (Fig. 2). This stems from two factors. 
Firstly, the inclusion of\J{l)e,n links in the parallel transporters increases the disorder in the 2-point 
correlation functions. In particular, the Z2 disorder in the phase of the U(l)e„, links for quarks leads 
to a Z2 electromagnetic flux string between constituent quarks in addition to the usual color flux 
string. What's more, since the explicit breaking of SU(2) center symmetry by quarks is suppressed 
by the introduction of fractional electric charge, the color links are less ordered than in standard 
2-color QCD. This is evident if one calculates meson masses using only the SU(2) links from our 
SU(2)xU(l)e,„ configurations. These masses, corresponding to electrically neutral sources, are 
consistent with those calculated on quenched SU(2) configurations at the same parameters. 

This begs the question: has the inclusion of fractional electric charge with respect to our U(l)em 
gauge action effectively removed the quarks from the functional integral? Such a conclusion would 
be reasonable if the Z2 \]{\)em phases for the quarks were truly random. But quarks loops induce 
a coupling between the SU(2) and U(l)e„, links. Consider the plaquette-plaquette and Polyakov 
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loop-Polyakov loop terms from our effective action Eq. (4. 1). These are minimized if the products 
are unity. That is, if the SU(2) and U(l)e,„ links order with respect to one another. This ordering is 
evident in the real part of the trace of the combined SU(2) xU(l)e,„ Polyakov loop that corresponds 
to a quark line 



-quark \ 



x) = cosPe/2(-?) • Re Pcoi{x) 



(4.2) 



which increases as fc is increased, i.e. as the quark mass decreases (see Fig. 3). 

We have explored the phase diagram with the ^(fc^) effective action Eq. (4.1), which allows 
us to cheaply check our intuition about the coupling of SU(2) and U(l)e,„ links without having to 
wony about the chiral limit for K. In Fig. (4) we show results for the SU(2) Polyakov loop with 
the U(l)e,M phases e'^''^^^/^ restricted to ±1, which amounts to jSe,„ = oo in the U(l)em gauge action. 
As K increases, the disorder-order transition of the SU(2) Polyakov loop moves to smaller values 
of j3co/ and sharpens dramatically. Note that in the combined limit K", jS^m — >■ °°, the plaquette- 
plaquette coupling in Eq. (4.1) forces the SU(2) plaquettes to take the values ±1. For very large K 
the transition line should therefore terminate with the first order bulk transition of Z2 gauge theory 
at j8fo/ ~ 0.44 [18]. As the SU(2) plaquette is driven to unity for lai^ge j3co/, U(l)e,„ plaquettes 
corresponding to quark loops, cos ^, receive an effective coupling that suppresses Z2 disorder in 
their phases. The relevant \]{\)e,„ Polyakov loop cosPg /2 remains disordered at small values of K, 
but transitions to unity for large ^coi and K (not shown). 

The coupling SU(2) and U(l)e„j links by quarks therefore leads to results distinct from the pure 
gauge theory, but the interesting physics lies at much larger values of K than for standard 2-color 
QCD. In particular, the mass scale has been radically shifted by the inclusion of fractional electric 
charge. A ratio of p and n masses of ~ 1.28, which is obtained for K = 0.178 at ^coi = 1.7 in 
standard 2-color QCD, requires K ~ 0.223 in our SU(2)xU(l)e,„ theory at j3e„, = 2. The chiral 
Umit is correspondingly pushed from Kc ~ 0.185 to ~ 0.241 (see Fig. 2). 
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Figure 3: Volume averaged SU(2)xU(l)(,,„ 
Polyakov loop corresponding to a static quark, 
Pquark{x) = cosPq/jC^) ' Re Pcoi{x), OH an 8^ X 4 lat- 
tice with jSfo/ = j5e,n = 2. SU(2) and U(l)„„ links 
order with respect to each other as the quark mass is 
decreased. 




Figure 4: SU(2) Polyakov loop on 16^ x 4 lattices, 
using an ff{K^) hopping effective action at jS^,,, = °°. 
The location of the transition is shifted with increas- 
ing K. 
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5. Discussion 

Given the fractional electric charge of quarks and the existence of a global symmetry that 
relates the centers of the color and electroweak gauge groups, it may be misleading to study non- 
perturbative phenomena in QCD alone. In our toy 2-color model, the coupling of Wilson quarks 
with half-integer electric charge relative to a compact \]{\)em has a dramatic effect on the color 
sector. The SU(2) Polyakov loop has an order-disorder transition that is consistent with the spon- 
taneous center symmetry breaking transition of the pure gauge theory at values of the hopping 
parameter where center symmetry is clearly explicitly broken in standard 2-color QCD. 

The model brings with it additional technical problems, however. Hybrid Monte Carlo trajec- 
tories must be supplemented with additional local updates to ensure ergodicity, and the mass scale 
is shifted to much larger values of K. The determination of a complete phase diagram will require 
considerable computational effort. As such, it is worthwhile to study fractional electric charge in 
even simpler models, such as our ^(k"^) hopping expansion model. Numerical simulations in an 
analogous gauge-Higgs theory are also ongoing [19]. 
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